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Abstract. The possibility of excitonic condensation in a recently proposed 
electrically biased double-layer graphene system is studied theoretically. The main 
emphasis is put on obtaining a reliable analytical estimate for the transition 
temperature into the excitonic state. As in a double-layer graphene system the total 
number of fermionic "flavors" is equal to iV = 8 due to two projections of spin, two 
valleys, and two layers, the large-iV approximation appears to be especially suitable 
for theoretical investigation of the system. On the other hand, the large number of 
flavors makes screening of the bare Coulomb interactions very efficient, which, together 
with the suppression of backscattering in graphene, leads to an extremely low energy 
of the excitonic condensation. It is shown that the effect of screening on the excitonic 
pairing is just as strong in the excitonic state as it is in the normal state. As a result, 
the value of the excitonic gap A is found to be in full agreement with the previously 
obtained estimate for the mean- field transition temperature Tc, the maximum possible 
value l^^^^^T^'^^ ~ 10~^eF (ef is the Fermi energy) of both being in ImK range 
for a perfectly clean system. This proves that the energy scale ~ 10~''ep really sets 
the upper bound for the transition temperature and invalidates the recently expressed 
conjecture about the high-temperature first-order transition into the excitonic state. 
These findings suggest that, unfortunately, the excitonic condensation in graphene 
double-layers can hardly be realized experimentally. 
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Figure 1. Excitonic condensate in a system of two spatially separated graphene layers. 
Electrons and holes in the layers are induced by applying the external gate voltage and 
the excitonic pairing is caused by the attractive interlayer Coulomb interaction. 

1. Introduction 

The idea of excitonic condensation in metallic systems was originally proposed [1] by 
Keldysh and Kopaev for semimetals with overlapping conduction and valence bands. 
They have shown that the attractive Coulomb interaction between electrons and holes 
leads to an instability towards formation of bound electron-hole pairs, analogous to 
the Cooper instability in superconductors. Somewhat later, Lozovik and Yudson 
suggested [2] that excitonic condensation could be realized in a double-layer system of 
spatially separated electrons and holes. If the layers are close enough to each other, the 
interlayer Coulomb interaction could still be appreciable, which would lead to coupling 
between electrons and holes. 

It took a while after the proposals of Refs. [2] until the technology able to fabricate 
electron-hole bilayers has been developed. So far, experimental efforts have been mainly 
concentrated on the GaAs/Ala;Gai_a;As double quantum well heterostructures and the 
evidence [3l HI [5l [6l [TJ [H [9], [10] for excitonic condensation in such systems based on the 
investigation of the Coulomb drag [HI HJ [5] and photoluminescence spectra [6l El [8], [9l [10] 
is building. Another exciting phenomenon in the coupled semiconductor bilayers is the 
quantum Hall ferromagnetism [TT| [12]. 

Besides semiconductor double quantum wells, what other materials and technology 
could be used for the search of the excitonic condensation? Graphene [I3l [TH [151 [161 
[T71 [T8l [T9] . an isolated atomic monolayer of carbon, could seem very attractive for 
this purpose. Indeed, quite recently several groups of authors [201 [2ll [221 [23] proposed 
the following graphene-based setup (Fig. [1]) as a candidate for the observation of the 
excitonic condensation. If one takes two graphene layers separated by an insulator, then 
electron doping in one layer and hole doping in the other can be obtained by applying 
external gate voltage. Analogously to the semiconductor quantum double wells, the 
interlayer Coulomb interaction would lead to coupling between electrons and holes. 

What attractive properties does such graphene system have? First, the relatively 
high values of the Fermi energy ep ~ 0.3 eV that can be achieved in graphene by 
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using electric gates [14j is an obvious advantage, since the condensation temperature is 
expected scale linearly with ep. Second, as any mismatch between the Fermi momenta 
of electrons and holes tends to suppress the condensate in the same way as the Zeeman 
splitting suppresses the conventional superconductivity, the nearly perfect matching of 
the electron and hole branches of the spectrum and the ability to fine-tune the carrier 
density become extremely important. 

So, in terms of engineering properties, the double-layer graphene system (Fig. [1]) 
could seem exceptionally attractive for realization of the excitonic condensation. 
However, before one proceeds with the fabrication of this system, an important 
theoretical question that should be answered is how high the temperature of transition 
to the excitonic state in such a setup could be. 

Some seriously conflicting predictions have been reported in the literature on 
this matter. Initially, an estimate for the Berezinski-Kostelitz-Thouless (BKT) 
transition [211 [25] temperature in a double-layer graphene system was made in 
Refs. [221 [23]. Studying the problem numerically, the authors came to conclusion that 
the BKT temperature T^j. in this system could be very large, 

T^T ~ O.le^. (1) 

With the Fermi energy ep ~ 0.3eV, this would correspond to room temperatures, 
T^j, ~ 300K, which would be very encouraging for experimentalists. 

Unfortunately, the estimate ([1]) for the transition temperature is too optimistic. 
As pointed out in Ref. [26] this result stemmed from the complete neglect of screening 
of the interlayer Coulomb interaction by the carriers in the layers in the analysis of 
Refs. [221 [23] . As the bare Coulomb interactions in graphene are not weak, it is not 
surprising that neglecting screening the transition temperature appeared to be only 
one order smaller than the Fermi energy ep- In Ref. [26] the mean field transition 
temperature Tc, at which the normal state becomes unstable towards electron-hole 
pair formation, was calculated. It was demonstrated that when screening is taken 
into account the highest possible value of the mean-field temperature turns out to be 
extremely small, 

T^max _ lQ-7^^_ (2) 

The reason for that is the large number = 8 of "flavors" in a double-layer graphene 
system, which makes screening very efficient, and the suppression of backscattering 
in graphene due to the chiral nature of Dirac quasiparticles. These factors combined 
make the system effectively weakly interacting, and the transition temperature depends 
exponentially on the coupling constant in the weak-coupling regime. Since the thermal 
fluctuations of the order parameter in two dimensions can only suppress the superfluid 
properties of the system, the BKT transition temperature T^t can be only smaller than 
the mean-field temperature ([2]). 

Somewhat later, in Ref. [27], the authors of Ref. [22] presented arguments to justify 
the neglect of screening in Refs. [221 123]. They argued that in the excitonic state 
screening had to be essentially suppressed due to the gap in the excitation spectrum. 
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This would mean that in the excitonic phase, at least at low enough temperatures, the 
interactions are not screened and the system is in the strong coupling regime of the bare 
Coulomb interactions. Therefore, the value of the gap A*(T = 0) at zero temperature, 
when the suppression of screening is maximal, had to be of the same order as the initial 
estimate ([I]) for the BKT temperature. 



A strong discrepancy between the small value of the mean-field temperature and 
the large value ([3]) of the gap in the excitonic state implies, according to the authors 
of Ref. [27], that the system should undergo a first order transition at temperature 
T ~ A*(T = 0), which would again correspond to room temperatures. 

The purpose of this paper is to review in more detail the arguments leading to 
Eqs. ([2]) and ([3]) in order to understand which of the estimates for the transition 
temperature is finally correct. We will show that, in terms of its effect on the electron- 
hole pairing, screening in the excitonic state is just as strong as it is in the normal 
state. This is because a wide range of scattering momenta of electron and hole in a 
pair contributes to the value of the excitonic gap, whereas screening is suppressed only 
in a much narrower range of small momenta. As a result, an accurate analysis of the 
excitonic state yields essentially the same estimate 



for the highest possible value of the gap as the one obtained for the mean field 
temperature [Eq. ([2])]. Therefore, the energy scale ~ 10~'^eF is the only one in the 
problem, no matter whether one considers the normal or excitonic state, and Eq. ([2]) is 
indeed the upper bound for the transition temperature into the excitonic state. This 
invalidates the argument of Ref. [27] about the high-temperature first-order transition. 
Considering that the excitonic condensate is sensitive to the impurity scattering, such 
a low transition temperature renders the observation of the excitonic condensation in a 
double-layer graphene system very improbable. 

The analysis of the problem presented in this paper is largely based on Ref. [28], 
in which a detailed theory of excitonic condensation in a single-layer graphene subject 
to the parallel magnetic field was developed. Although here we mainly concentrate 
on obtaining an analytical estimate for the transition temperature, many of the fine 
features of the excitonic state studied in Ref. [28] could be carried over to the double- 
layer graphene system. 

2. Large- approach to a double-layer graphene system 

In this section, we will outline the main idea of the theoretical method, the large- A^ 
approximation, which will be used to obtain analytical estimates ([2]) and (jl]) for the 
transition temperature and the excitonic gap. For a double layer graphene system the 
large-A^ approximation appears to be particularly reliable and is expected to provide 
good quantitative predictions. 



A*(r = 0) ~ O.Uf. 



(3) 



^max^^ = 0) ~ lO'^ep 



(4) 
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Figure 2. The idea of the large- approximation. In each order in the bare Coulomb 
interaction (wavy lines), the leading in A'' ^ 1 diagrams for the effective interaction 
( "photon propagator" , double wavy line) are those that contain the maximum number 
of fermionic loops. So, in the third order [diagrams (3.1)-(3.3)], the diagrams (3.2) 
and (3.3) are smaller than (3.1) in 1/N <C 1 and can be neglected. Performing the 
same analysis for higher orders, one arrives at the series analogous to that of the RPA. 
The summation of this series can be reduced to solving an equation for the effective 
interaction (last line). 

The bare strength = j[ev) (e is the electron charge, e is the dielectric constant 
of the insulating medium surrounding graphene, v ^ lO^cm/s is the velocity of the Dirac 
spectrum, and we will use the units, in which the Planck's constant h = 1 throughout 
the paper) of the Coulomb interactions in graphene is not that small. For graphene in 
vacuum {e = 1) 2.16 and for Si02 as an insulator typically 1 or somewhat 
smaller. These values of are not exceptionally large, but cannot be considered as 
small either. Therefore, strictly speaking, the weak-couphng treatment of the Dirac 
fermions in graphene for the realistic values of Tg is not applicable. In the diagrammatic 
approach, for > 1 one would have to sum up all diagrams in each order in the 
bare interaction. Clearly, this task cannot be carried out analytically. If, however, the 
number N of independent fermionic species, or "flavors" is large, one has an additional 
expansion parameter, which can efficiently be exploited. The flavors can correspond to, 
e.g., different projections of the real spin S of particles, in which case would be the 
total number of possible spin projections, N — 2S + 1, or to some other effective discrete 
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degrees of freedom, which stem from the properties of the particles spectrum. 

So, if the number of flavors N is indeed large, one can notice (Fig. [2]) that after 
summing the contributions of all flavors ( "performing the trace over the flavor indices" ) 
each closed fermionic loop consisting of the fermionic Green's functions gives the factor 
A^. Therefore, in a given order in the bare interaction, the leading in ^ 1 contribution 
to the effective interaction ("photon propagator") comes from the diagrams with the 
maximum numbers of electron loops. These are the diagrams that contain only the 
polarization bubbles consisting of two Green's functions. The series of such diagrams 
can already be summed up as it is formally identical to the one of the well-known 
Random Phase Approximation (RPA). 

Physically, the leading in A^ series describes the effect of linear screening of 
interactions by the carriers. Of course, screening is to some extent present in any system 
of free charged particles. However, in the limit of large A^ screening is particularly strong, 
since all N species participate in the screening of the interactions between fermions 
of each particular species. Screening makes the system effectively weakly interacting 
reducing the coupling constant from to 1/A^: 1/N. This allows one to study 

the system analytically, despite the fact that the bare Coulomb interactions may be not 
weak (r, > 1). 

How large is the number of flavors A^ in a double-layer graphene system (Fig. [1])? 
In a single layer of graphene the number of species is equal to 

N, = NsN, = 4 

due to A^s = 2 possible projections of spin and = 2 valleys. The large- A^ approach to 
a single-layer graphene was used in a number of works [281 ESI EOl EH [321 [331 [311 [33 [36] 
before. For a double-layer system one has an additional "which-layer" degree of freedom, 
since each carrier can belong to either one of the layers. This twofold {Ni = 2) 
degeneracy appears formally in calculations when summing the polarization bubbles 
of both layers. Therefore, in a double-layer graphene system (Fig. [1]) the total number 
of flavors is equal to 

N = NsN,Ni = 2^ = 8. (5) 

For this value it is quite reasonable to expect the large- A^ approximation to give good 
quantitative predictions. The accuracy of the obtained results will be estimated below 
in Sec. [H 

So, with N = 8 the large-A^ approximation seems to be particularly suitable for 
the double-layer graphene system. On the one hand, this is an attractive situation 
for theorists, since one is able to perform analytical calculations and make reliable 
quantitative predictions for a system of non-weakly interacting electrons. On the other 
hand, since the system is effectively in the weak coupling regime, one could expect the 
transition temperature to be quite small, which certainly is a pessimistic prospective for 
experimentalists trying to observe excitonic condensation in such a system. 

Before we proceed with calculations, we would like to emphasize one important 
point. Although the diagrammatic series in RPA and large- A^ approximation are 
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formally identical, there is an important difference in the validity of the two approaches. 
The RPA works, when the characteristic energies e and momenta q (set, e.g., by 
temperature) of the excitations contributing to the studied quantity are small compared 
to the Fermi scale: e <^ ep and q <^ pp. The reason for that is the possibility to describe 
the low-energy excitations in terms of the charge and spin densities, quadratic in the 
fermionic operators. One the other hand, the large-A^ approximation is applicable as 
long as ^ 1 and does not require the smallness of the relevant energy scales compared 
to the Fermi energy. It works even when the Fermi surface is absent in the problem, e.g., 
in the vicinity of the Dirac point in graphene [28l[29l[30l[3ll|32l[33l[3l[35l[36]. Actually, 
as we will see in the next section, the excitonic gap and the transition temperature 
are determined by all transfer momenta q in the range < g < 2pF, where pp is the 
Fermi momentum. Therefore, strictly speaking, the RPA does not apply to the problem 
when the bare interactions are not weak. It does, however, apply to the case of weak 
bare interactions and we will come back to this point, when we discuss the transition 
temperature in more detail in Sec. HI 



3. Gap equation 

In this section, we present the self-consistency equation for the excitonic gap, which will 
be used in the next two sections to obtain the mean-field temperature and the gap in 
the excitonic state. 

To study the excitonic condensation in a double-layer graphene system (Fig. [1]) we 
will use the following electron Hamiltonian, 

i=e,h i,j=e,h 

Here, the indices i and j numerate electron (e) and hole (h) layers. The terms 

= J dr0l(r)(t;Tp-/x.)0.(r) (7) 

describe the free Dirac particles in graphene layers, where 0i(r) are the Dirac spinor 
field operators, r is a two-dimensional radius vector in graphene layers, p = — iVr is the 
momentum operator, and r = (r^, Ty) is a vector consisting of the Pauli matrices and 
Ty in the sublattice space of graphene honeycomb lattice. The chemical potentials fie,h 
of electrons and holes have opposite signs and will be assumed equal in absolute value, 

/ie = —fih = 

This corresponds to the most favorable situation for the excitonic condensation, as any 
mismatch /Xg— |yU/i| acts as an effective Zeeman field and tends to suppress the condensate. 
Further, the terms 

Ht = \j dr dr'0l(r)0j(r')r.,(r - r')0,(r')0.(r) (8) 
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in Eq. describe the intralayer {i = j) and interlayer {i ^ j) Coulomb interactions 
between the Dirac particles with Vo(r) = Vee(r) = V/ift(r) and Vi(r) = Veh{r) = Vhe{r) 
given by the Coulomb potential 

= ^, V^i(r) = -J==r,^ (9) 
|r| ^|r|2 + c/2 

In Eq. ([9]), is an effective electron charge screened by the insulator embedding 
graphene sheets, el = e'^/e, and d is the distance between the layers. 
The order parameter of the excitonic condensate is defined as 

A(n) = J dre-i^'^"(^-'-Vi(r - r')(0e(r)0i(r')). (10) 

Here, n is the unit vector representing the direction of the electron momentum on the 
Fermi surface. In the sublattice space A(n) has the following matrix structure, 

A(n) = AP(n) + A'P(-n), (11) 

where V{n) = (ir^ + Tt)/2 and t = {—ny,nx) is the vector orthogonal to n = {rix^ny) 
in the xy plane. The constant A is the excitonic gap of the spectrum, = a/AM-~^, 
^ = v{p — pf). The physical quantities of the condensate are determined by A, whereas 
the other constant in Eq. ffTTl) . A', which is related to A, does not enter any of them. 
The properties of the excitonic gap A will be of our main interest from now on. 

Within the large-A^ approach, the theory is analogous to the conventional weak- 
coupling Bardeen-Cooper-Schriffer (BCS) theory with the interaction given by the 
screened interlayer Coulomb potential. The equation for the excitonic gap A has the 
conventional BCS-like form [28j 

n=—oo " '^t II ^ 

In Eq. ( fT2l) . T is the temperature, = TTT{2n + 1) are the fermionic Matsubara 
frequencies and the summation should be carried out over all integer n. The integration 
over ^, which diverges logarithmically if the summation over n is performed first, should 
be cut by the Fermi energy ±6^. Further, the dimensionless coupling constant of the 
screened interlayer Coulomb interactions is given by the expression 

, ( . e\ i + cosg 

Xa = i^ 7;-»^a SpiT'Sm- , (13) 







27r V 2 y 2 

where u = ep/{2nv'^) is the density of states in each graphene layer per one projection 
of spin and one valley. In Eq. f[T^ . V^{q) is the Fourier component of the static 
(at frequency u = 0) screened interlayer Coulomb interaction potential. The static 
limit for the coupling potential is valid with the logarithmic accuracy. The potential 
V/^{q) depends on the value of the gap A, since screening is affected by the presence 
of the excitonic gap. We will discuss the latter point in detail and present the explicit 
expression for VA{q) in the excitonic phase in Sec. [5l 

Two important points should be emphasized about Eqs. (fT2|) and (fT3|l . First, as 
follows from Eq. (|T3l) . the coupling strength Aa is determined by all transfer momenta 
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q = 2pp sm{9/2) in the range < g < 2pp, or equivalently, by all scattering angles in 
the range < 6 < 2tt. This integration over the scattering angles is standard and it 
also appears in the conventional BCS theory. This fact is, however, crucial for obtaining 
the correct estimate for the gap A in the excitonic state. It will be key to disproving 
the statement made in Ref. [27] that, in terms of its effect on the magnitude of the gap, 
screening is suppressed in the excitonic state. 

The second point is the presence of the factor (1 + cos 9)/ 2 in Eq. f|T3|) . It stems 
from the chiral nature of Dirac fermions in graphene and leads to the suppression of 
backscattering. This feature is specific to graphene and is absent in the system with 
conventional electronic spectrum. So, for any point-like interaction, Vi^{q) = Uo = const, 
the coupling strength Aa = vUq/2 is actually two times smaller than the conventionally 
defined dimensionless coupling constant vUq. As the transition temperature following 
from Eq. (fT2l) depends on Aa exponentially (see next section), this reduction by the 
factor 2 leads to an additional substantial suppression of the transition temperature 
compared to the double layer system with the conventional electronic spectrum. 

Equations f|T2l) and f|T3|) supplemented by the proper form of the interlayer Coulomb 
potential Va(q') allow one to obtain an analytical estimate for the transition temperature 
and the excitonic gap. In the next section we will calculate the mean field transition 
temperature. The results of the next section reproduce those of Ref. [26] . 

4. Mean-field transition temperature 

In order to find the mean field transition temperature Tc, one needs to linearize the gap 
equation ( fT2l) in A and determine the temperature at which a nonzero solution to it 
appears. Calculating the integral over ^ and the sum over n in the linearized Eq. (fT2l) . 
one arrives at the BCS-like expression 



In Eq. dHj), the screened Coulomb potential Va=o(q') [see Eq. ( IT3i) for Aa=o] should be 
determined in the normal state, when A = 0. 

Calculating the potential Va=o(?) in the large- approximation is formally 
analogous to the RPA. One important difference, however, is that since the relevant 
transfer momenta q = 2pp sin ^/2 belong to the range < g < 2pp [Eq. f|T3|) ]. one has to 
use the exact expression for the polarization operator 11 a (t^, g), and not its approximate 
form at g -C pp. However, the static polarization operator nA=o('?) = nA=o(i^ = 0,g) 
in the normal state in graphene [Ml [37] does not depend on momentum q at all in this 
range and equals 



The expressions for the bare intralayer and interlayer Coulomb interaction potentials 
[Eq.(l9])] in the Fourier representation read 



Tc ^ exp(-l/AA=o)eF- 



(14) 



nA=o(g) = NsN^v, q < 2pp. 



(15) 



Vo{q) 



27161 exp(— grf) 



(16) 
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respectively. Using Eqs. (fTB]) and (ITBl) and performing analogous to the RPA calculations 
one obtains for the screened interlayer Coulomb potential in the normal state 

T/ f \ 27relexp{-qd) 

g + zx + — exp(— zgajj/g 



In Eq. (d 

X = 27rA^,A/'^e^z/ (18) 

is the Debye screening wavevector (inverse screening length) in each layer. 

Substituting Eq. f|T71) into Eq. f|T3l) . one obtains the mean field transition 
temperature Tc from Eq. ( |T4|) . Let us estimate the highest possible value of Tc. 
According to Eq. (ITTl) the potential Va=o(q') is maximal, when the interlayer spacing 
d satisfies the condition 

xrf < 1. (19) 

The other required condition, ppd -C 1, is automatically satisfied for > 1 in the 
limit (flQll . since x = rgNgNypp = '^TsPf- In reality, the limit ( |T9|) is hardly achievable. 
Indeed, on the one hand, one would like to have a large Fermi energy ep to have a higher 
cutoff energy in Eq. (fT4l) . but on the other hand, higher ep would require very smaller 
interlayer distance d. For ep = 0.3eV and = 1, the condition (fT9l) corresponds to 
d -C 0.5 nm. 

So, in the limit ( |T9l) of "zero" interlayer distance (i, the interlayer potential ( |T71) is 
maximal and equals 

Equation (120|) has the same form as the screened Coulomb potential of a single layer. 
Note, however, the factor 2 in the term 2x in Eq. ( I20l) . It originates from summing 
the polarization operators of both layers and reflects the above discussed "which-layer" 
degeneracy {Ni = 2), while the Debye wavevector x [Eq. ( |T8l) ] contains the spin (A^^ = 2) 
and valley (A^„ = 2) degeneracies of each layer. One may say that screening in a double- 
layer system is twice as efficient as in a single layer, since the carriers of both layers 
participate in it. The potential V^Q(g) is a decreasing function of q with the maximum 
at g = 0. Using Eq. f|T8l) we obtain 



uVI^Uq = 0) = ^ = ^- (21) 

We see that the maximum fl2Tl) of the dimensionless coupling constant is given by the 
inverse number of flavors and is, therefore, small for large A^. Inserting Eq. fl^Tl) into 
Eq. (fT3|) . one obtains the maximum possible value of the coupling constant Aa=o in the 
normal state: 

XT^ = — = L_ (22) 

As discussed above, the coupling constant X^To^ which determines the transition 
temperature, is twice as small as the one given by Eq. (12T|) due to the suppression 
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of backscattering. Inserting Eq. fl22p into Eq. (1141) . one obtains the maximum possible 
value of the mean-field transition temperature, 

T,'"^" ^ exp{-2N)eF = exp(-16)e^ ^ 10-^- (23) 

We see that, although the value X^To [Eq- ( |22l) ] is not particularly small, the mean-field 
temperature appears to be extremely small due to its exponential dependence on A^^q. 
With ep ~ 0.3eV, Eq. ([23]) corresponds to T^'"^^ ~ ImK. 

Let us here emphasize the difference between the large- approximation and the 
RPA in connection to the considered problem. Basically, Eqs. f|T2l) and f|T3l) of the 
weak-coupling BCS-like theory are valid whenever the dimensionless coupling constant 
Aa of the screened interlayer Coulomb interactions is small. In case of a sufficiently large 
number of flavors the coupling constant A^t^o [Eq- ( |22l) ] is small in <C 1 and the 
applicability of the weak-couphng theory [Eqs. f|T2l) and f|T3|) ] is justified for arbitrary 
strength of the bare interactions. However, if the number of flavors were not large (e.g., 
if A^<j = A^^ = 1), Eqs. ( fT2|) and ( |T3l) would still be valid in the weak coupling regime 
of the bare interlayer interaction Vi(g), i.e., when ^ 1 or ppd ^ 1, see Eq. ( IT6|) . 
Indeed, in the integral over q = 2pFsm{9/2) in Eq. (|T3l) . z/Vi(g ~ pp) ~ <^ 1 in 
the major range of integration q ~ pp, and for these q the bare interaction potential is 
weak. Only for q ^ pp the bare coupling vVi{q -C pp) ^ 1 is not small. But for small 
q <^Pp the RPA is valid and the effective interaction is described by Eq. (fTTl) . in which 
the number of flavors N does not need to be large. Actually, the gap equations (fT2|) . 
(|T3l) . and ( ITTl) were obtained in Refs. [20l [21] in limit of weak interlayer interactions 
(r^ <^ 1 or p^d <^ 1) within the RPA approximation before Eq. f|T^ for the mean-field 
temperature was derived in Ref. [26] using the large- iV approximation. An important 
conclusion of Ref. [26] is that due to the reasonably large numerical value = 8, these 
equations are valid and provide reliable quantitative prediction for not only for weak 
interlayer interactions, but also in the realistic case of moderate to strong interactions, 
when > 1. 

Let us estimate the accuracy of the result ( l23i) . This equation is obtained by taking 
into account the leading in A^ diagrams (Fig. [2]) for the effective interlayer interaction. 
The diagrams that are neglected [such as diagrams (3.2) and (3.3) compared to (3.1)] 
are smaller in 1/A^ than the leading contribution. Therefore, one could represent the 
exact coupling constant A^^^q in the form of an asymptotic series in A^ with Eq. fl22l) 
giving the leading term of the series. This way, one can write down the inverse coupling 
constant as 

l/AS5) = 2Ar-C»(l), (24) 

where 0{1) > represents all subleading in 1/A^ terms and is of order unity, 0{1) ~ 1. 

Inserting Eq. (|2^ into the exponential function in Eq. (|T^ one can conclude that 
in the best case scenario the transition temperature T™*^^ could be about one order 
(^gC>(i)) greater than the one given by Eq. ( l23i) . With the proper (yet unknown within the 
logarithmic accuracy) numerical prefactor in Eq. (123|) one could optimistically hope that 
jimax QQ^y reach lOOmK for ep ~ 0.3eV. Unfortunately, what concerns the possibility to 
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observe the effect experimentally, this would not radically change the situation. First, as 
it was mentioned, the maximum fl22l) of the interaction strength is virtually unachievable 
in the experiment, since this requires the condition ( |T9|) to be satisfied. Second, as it will 
be discussed in Sec. [6l inevitable boundary scattering at the edges of graphene sheets 
will destroy any condensate with Tc < IK even in a system free of any bulk disorder. 

5. Excitonic gap and screening in the excitonic state 

In addition to the mean field temperature Tc obtained in the previous section, one 
can estimate the Berezinski-Kosterlitz-Thouless (BKT) transition temperature. The 
BKT transition p4l [25] occurs in two-dimensional superfluid systems, in which thermal 
fluctuations destroy the long-range order of the order parameter. The BKT temperature 
is defined as the temperature at which creation of topological defects (most commonly, 
vortices) in the order parameter becomes thermodynamically favorable. Mathematically, 
this condition can be expressed as: 



where p(T) is the so-called stiffness of the order parameter, which determines the energy 
of the topological defect. The factor Nkt = 1 for vortices. As shown in Ref. |28j, 
however, in graphene creation of halfvortices can be more favorable, in which case 
Nkt = 4 and Tkt would be somewhat lower than for Nkt = 1- Within the mean- 
field theory, the stiffness has the form [28] 



where A(T) is the excitonic gap, which should be obtained by solving the gap equation 
(IT^ . The stiffness (12^ takes a universal value p(T = 0) = ep/{4:7i) at low temperatures, 
when T <ti A(T), decreases with increasing temperature, and vanishes when T 
approaches the mean-field temperature (fn|) . T — > T^, where the gap turns to zero, 
A(Tc) = 0. Therefore, the BKT temperature Tkt, which is a solution of Eq. (125|1 . 
cannot exceed the mean-field temperature T^, 



This is illustrated graphically in Fig. [3l The result (1271) is also intuitively clear, since 
one would naturally expect thermal fluctuations of the order parameter to suppress, 
rather than enhance, the superfluidity of the system. 

One can conclude from Eq. fl27|) that the estimate ([2]), obtained within the mean- 
field approach, sets the upper bound for the temperature of any realistic [28] phase 
transition into the excitonic state in a double-layer graphene system. Nevertheless, 
as mentioned in the introduction, the authors of Ref. [27] came to quite a different 
conclusion. On the one hand, they agreed that the mean-field temperature Tc indeed 

I The authors of Ref. [27] actuaUy provided an estimate Tc ~ 6-10~^eF for the mean-field temperature, 
which is more than two orders greater the result ([2|) for the highest possible value of Tc. It is hard to 



p{Tkt) — NktTkt, 



(25) 




(26) 



Tkt < Tc. 



(27) 
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Tkt Tc 

Figure 3. Graphical solution of Eq. (|25p for the BKT transition temperature 
Tkt- The stiffness p{T) [Eq. ([^ J of the order parameter takes a universal value 
p{T = 0) = eF/(47r) at low temperatures (when T ^ A(r)) decreases with increasing 
temperature T, and vanishes at the mean-field transition temperature T = Tc, where 
the excitonic gap A(T — Tc) = turns to zero. Therefore, the BKT temperature Tkt 
cannot exceed the mean- field temperature Tc, Eq. ([27]) . 



appears to be significantly smaller than their initial estimate T^t [-^1- ©] BKT 
temperature, provided screening is taken into account. On the other hand, they argued 
that in the excitonic phase the gap in the excitation spectrum suppresses screening, 
which would mean that the gap A*(T) itself should be of the same order as the 
estimate ([I]) obtained neglecting screening, see Eq. ([3]). Such a strong discrepancy 
between [Eq. ^] and A*(T = 0) [Eq. ©] indicates, according to Ref. [27j, that the 
system should undergo a first-order transition into the excitonic state at temperature 
T ~ A*(T = 0), which would again correspond to high temperatures. 

We will now show that, in fact, the other, less exciting, scenario is realized. We 
will demonstrate that the suppression of screening by the gap in the excitonic state, 
although it does take place, has a negligible effect on the value of the gap, when the 
latter is determined self-consistently from the gap equations ( fT2l) and ( fT3l) . As a result, 
the system is in the weak-coupling regime of the screened interactions, described by 
the large-iV approximation, at all temperatures below and the value of the gap 
A(T) < 10~^eF in the excitonic state appears to be in full agreement with the value ([2]) 
of the mean-field temperature. 

To prove this point, one needs the expression for the screened interlayer Coulomb 
potential V/^{q) in the excitonic phase, where A 7^ 0. Performing the calculations 
analogous to those of RPA, in the large- approximation one obtains 

y ^ V-jq) .281 

^^^^ l + 2n+(g)W(g) l + 2U-^{q)V4qy ^ ' 

Here 

1 27re^ 
V±{q) = -[l±exp{-qd)]-f (29) 
z q 

comment on the possible origin of this discrepancy, since no analytical formula for Tc was presented in 
Ref. [H]. 
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f+oo „2 
1-nT J2 J 



are the half-sum and half-difference of the intralayer Vo{q) and interlayer Vi{q) bare 
Coulomb potentials [Eq. f|T6|) ]. The polarization operators 

ni(g) = n°,(g)±ni(g) (30) 

in Eq. ( l28l) are the sum and difference of the intralayer Il%{q) and interlayer Il]^{q) 
polarization operators. The operators Il%{q) and Il]^{q) describe the density response 
in one layer to the density perturbation in the same and other layer, respectively, and 
can be written down in the Matsubara representation [38] as 

n^'(9) = I'^^dr J dre-''^(A(r,r)p,(0,0)) (31) 

with the same layer indices {i = j) for Il%{q) and different ones {i ^ j) for 11^ (g). In 
Eq. ( 13T1) . Pi(r) = 0j(r)0j(r) are the density operators and r is the Matsubara time. 

The polarization operators (13T|) can be calculated using the Green's functions for 
the excitonic state. Doing so, one obtains that in the excitonic state 11^(5') does not 
depend on the gap at all and is equal to its normal state value 

n^(g) = N,N,u. (32) 

On the hand, the operator n^(q') is affected by the gap and the expression for it reads 

r^^de_ A2 

2^ [{el + A2) + {vqcosey/A]y/^^TX^ 

Here, as in Eq. ( fT2l) . £„ = 7rT{2n + 1) are the fermionic Matsubara frequencies and the 
summation is done over all integer n. 

Aiming to find the highest possible value of the gap, let us again consider the limit 
(fT9|) of "zero" interlayer distance, when the potential V^{q) [Eq. (!28|) ] is maximal. In 
this limit the second term in Eq. (l28il can be neglected, V+{q) = Vo(g) [Eqs. (fT6ll and 
( 129|) ]. and the interlayer interaction potential takes the form 

1 + 2nx(g)K)(g) 

Again, as in Eq. (1201) . the factor 2 in the denominator in Eq. (!34|) originates from 
summing the polarization operators of both layers. 

Let us discuss the properties of the polarization operator fj33l) and the corresponding 
properties of the potential ( IMI) . In the normal state (A = 0), the second term in Eq. ( l33l) 
is zero, n^^Q(g) is given by Eq. (fT5|) and one recovers the expression (!20|) for V^^{q) 
from Eq. The second term in Eq. (!33|) exists in the excitonic phase only, when 

A 7^ 0. It is a decreasing function of q and approaches zero for q ^ A/v exceeding the 
condensate momentum A/v. So, for momenta q ^ A/v the polarization operator 

n+(g > A/v) ^ NsKv (35) 

is given by its normal state value ( |T5l) at any temperature, including T = 0, and the 
interlayer interaction potential is completely screened, 

^^max^^ > A/v) ^ ^. (36) 
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Figure 4. Dependence of the polarization operator n^(g) [Eq. (1551) . 1*3^ graph] 
and the maximum of the screened interlayer potential V^^^{q) [Eq. (IM)) . right graph] 
on momentum q in the excitonic state. The case of A/ep = 10~^, = 1, and 
zero temperature is shown. The operator nj(q) is suppressed for q < A/v below 
the condensate momentum A/u, whereas for greater q 3> A/w it approaches the 
normal-state value IlA=o{q) = vNgNy [Eq. ([T5|), dashed hne]. Consequently, the 
potential vV^^^{q) is unscreened and exceeds its small normal-state value vV^^{q) < 
1/8 [Eqs. (pOI and (pij) ] only for q < A/v. As the excitonic gap A is determined by 
all momenta in the range < q < 2pF, see Eqs. (IT^ and (fT5|) . the suppression of 
screening at g < A/u gives only a minor correction to the coupling strength Aa , which 
is small in the ratio A/ep [Eq- ([55]) ]. 



Therefore, an important conclusion that can be drawn from Eqs. (1351) and fl36|) is that 
for momenta q exceeding the condensate momentum A/v screening is just as strong in 
the excitonic state as it is in the normal state. 

At the same time, the situation is different for small momenta. The polarization 
operator n^(g) is indeed suppressed at g < A/t>, when the second term in Eq. ([2SD 
contributes considerably. The smaller q and T are, the stronger the suppression of 
nj(g) is. In fact, at zero temperature T = and zero momentum q = the second 
term in Eq. (133!) cancels the first term exactly and one gets 

ni(g = 0)|T/A=o = 0. (37) 

This means that at small momenta and low temperature screening is indeed suppressed 
and the interlayer Coulomb potential is unscreened, 

27re2 

Vr''{q<A/v)\T«A^^. (38) 

Q 

Equation fl37|) has a simple physical interpretation: at zero momentum and temperature 
all electrons and holes are paired into dipoles and there are no free carriers (excitations) 
that could screen the interaction. 

The dependence of the polarization operator n^(g) [Eq. (l33ll ] and of the interaction 
potential V^^^{q) [Eq. flMl) ] on q, illustrating the above properties, is plotted in Fig. HI 

It is essentially the property (|37|) of the polarization operator fl33|) that the authors 
of Ref. [22] used in Ref. [27] to argue that, in terms of its effect on the value of the 
gap A, screening should be suppressed. The expression (133|) at g = can be rewritten 
in the form of Eq. (4) of Ref. [27] using the conventional method [38] of replacing 
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the sum over the Matsubara frequencies e„ by the integral over a continuous energy 
variable. According to Ref. [27], the property fl37j) justifies the neglect of screening, 
when determining the value of the gap in the condensed phase. 

It is easy to understand now from Eqs. (fT2|) . ffT3l) . (l33l) and and the discussed 
properties of the polarization operator (1331) that the suppression of screening at small 
momenta q < A/v has a very little, in fact, negligible effect on the value of the gap. 
Indeed, as follows from Eq. f|T3|) and was pointed out earlier in Sec. [3|, the coupling 
strength A a and, consequently, the gap A [Eq. ( |T2l) ]. is determined by all scattering 
momenta q = 2pp sin6'/2 in the range < g < 2pF. As we have just discussed, screening 
is suppressed by the gap in the range < g < A/t> only [Eqs. (1371) and (1381) ]. whereas 
in the range A/v ^ q <2pp screening is just as strong as in the normal state [Eqs. fl35l) 
and ( l36l) . Fig. H]. Therefore, the difference 

5Aa = Aa — Aa=o 

between the pairing constant Aa [Eq. ( fT3l) ] in the condensed phase and its value Aa=o 
in the normal state comes from the range of momenta < g < A/t; only and is small 
in the ratio of the gap A to the Fermi energy ep = vpp, 

5Aa/Aa=o ~ A/ci.. (39) 

Equation ( |39l) shows unambiguously that the presence of the excitonic gap does 
not influence the value of the coupling constant as long as A ^ e^?. In other words, the 
statements of Ref. [27j cannot be correct unless the excitonic gap becomes of the order 
of the Fermi energy. 

Let us now flnd the maximum value A™'^^ [in the limit (fT9|) of "zero" interlayer 
distance] of the gap Aq = A(T = 0) at zero temperature, when the suppression of 
screening by the gap is greatest. From Eq. (|T2|) . one obtains that Aq satisfles the 
equation 

Ao ^exp(-l/AAo)eF. (40) 
Unlike Eq. (fT4l) . which provides the explicit expression for the mean fleld transition 
temperature, one has to solve Eq. (HOjl to obtain Aq, since Aaq depends on Aq. However, 
considering the extremely small numerical value of T^^^/ep [Eq. (1231) ] we conclude 
immediately that 

A^" ^ exp{~l/X2To) = exp(-16) ep ^ 10"^^, (41) 

for which the suppression of screening by the gap is neglected completely, is the solution 
of Eq. (HOl) . Indeed, inserting the value (HTl) into Eq. ( |T3i) and calculating the integrals 
numerically, we obtain that SX'^^^/Xa=o ~ 10"^ [Eq. (!39|) ]. and therefore A^^^ is 
determined by Eq. ( 1411) with the relative precision ~ 2N6X^^^/Xa=o ~ 10~^. 

Besides Eq. pTl) . there are no other solutions to Eq. fj40|) in the range Aq < ep, 
which can be checked numerically. One could use the highest hypothetical value ([3]) of 
the gap as the zeroth approximation, insert it into the right-hand side of Eq. ( l40l) . and 
obtain that already on the second iteration the solution collapses to the value (14T]) with 
high precision. 
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Figure 5. Sensitivity of the excitonic condensate to the impurity scattering. 
Impurities with the size of potential smaUer than the interlayer distance scatter 
electrons and holes not identically, thereby breaking electron-hole pairs and suppressing 
the condensate. 

We conclude that, contrary of the arguments of Ref. [27], in terms of its effect on 
the excitonic pairing, screening is as strong in the excitonic state as it is in the normal 
one. Partial suppression of screening at small momenta q < A/v has a negligible effect 
on the value of the gap when the latter is determined self-consistently from the gap 
equation. This clearly follows from the fact that all the scattering momenta q in the 
range q < 2pF contribute to the value of the gap, and not only small ones. Therefore, 
the normal-state expression f|T7|) for the screened intrelayer Coulomb potential can be 
used with great precision to study the excitonic phase at any temperature T <Tc below 
the mean-transition temperature ([2]). As a result, the maximum possible value of 
the zero-temperature gap (14T1) appears to be in full agreement with the result ([2]) for 
the mean field transition temperature. The fact that exp(— 16)e^ ~ 10~'^eF is the only 
energy scale of the excitonic pairing in the double-layer graphene system proves that 
Eq. ([2]) really sets the upper bound for the transition temperature and invalidates the 
conjecture of Ref. [27] about the high-temperature first order transition. 

We would like to stress that all the above discussion about how screening is affected 
by the gap in the excitonic state was not actually in any way specific to graphene. The 
very same analysis could have been carried out for conventional systems [H [2] with 
quadratic electron spectrum. What concerns graphene, the point that screening is only 
negligibly affected by the gap has been made earlier in Ref. [28|, where the dependence 
of the polarization operator on A in the excitonic state (its suppression in the minor 
range of momenta q < A/v) was neglected in the calculations for the reasons discussed 
above. 

6. Influence of disorder 

Although the estimate ([2]) predicts a really small value of the transition temperature, 
one could still hope to observe excitonic effects in a ImK range. Unfortunately, in a 
realistic system even at these low temperatures the excitonic pairing could hardly be 
realized. The reason for this is that these estimates were made for a system that is free of 
any kind of disorder, while the excitonic condensate in double-layer systems is sensitive 
to the impurity scattering [391 Ej- Indeed, since the bound electron and hole carry the 
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same momentum p, any scattering process that changes the momentum of electron and 
hole not identically, i.e., p ^ Pe for electron and p Ph for hole, so that Pe 7^ p^, 
breaks the electron-hole pair (see Fig. [5]). This is the case for impurities with the range 
of the scattering potential less than the interlayer distance d, since the potential of such 
impurities differs in the two layers. Therefore, any sufficiently strong-range impurities 
destroy the excitonic condensate. 

Originally, the effect of the impurity scattering on the excitonic condensate 
was studied analytically for conventional systems in Refs. [39l [2]. The theory is 
formally analogous to the Abrikosov-Gorkov's theory [IQ] for magnetic impurities in 
superconductors. For a double layer graphene system, analogous calculations were 
performed in Ref. [H]. The main result is that sufficiently short-range impurities with 
the scattering time r destroy the excitonic condensate completely as soon as 

rTc < 1, (42) 

where is the transition temperature of the ideally clean system. Equivalently, for 
the condensate to exist, electron momentum has to be conserved at the scale of the 
correlation length v/Tc. How weak impurity scattering in a realistic system can be? 
Even if the system is free of bulk disorder, the mean free path is limited by the sample 
size due to the boundary scattering at the sample edges. With the typical size of 
graphene devices ~ l/im this corresponds to the scattering rate 1/r ~ IK. Therefore, 
we can conclude that any condensate with Tc in a clean system below IK would be 
completely destroyed in a realistic system. 

7. Comparison with conventional double-layer systems 

In terms of engineering properties, the double-layer graphene system could seem as 
an exceptional candidate for the realization of the excitonic condensation. However, 
as we have seen in the previous sections, the peculiar properties of graphene spectrum 
(additional valley degrees of freedom and chirality) appear to play a negative role for this 
effect. To illustrate this point, let us estimate for comparison the highest possible value 
of the transition temperature in a double-layer system with the conventional quadratic 
electron spectrum, such as a semiconductor double quantum well, using the large- 
approximation. 

How should the equations of Sec. [3] be modified in this case? First, since the valley 
degeneracy = 2 is absent, the total number of flavors is = A^gA"; = 4 due to the 
two [Ng = 2) projections of spin and two layers (A^; = 2), Therefore, the maximum of 
the dimensionless coupling constant equals 

uV^:Uq = 0) = ^ = i (43) 

Second, the factor (1 + cos0)/2, which stems from the chiral nature of Dirac 
quasiparticles and results in the suppression of backscattering, is absent in Eq. (fT3|) . 
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Therefore, the maximum couphng strength A^^^q that determines the gap according to 
Eq. ([12]) is just equal to Eq. (gS]), 



Ar=o = ^ = i- (44) 

As we see, the pairing strength A^^q in a double-layer with conventional spectrum is four 
times greater than the one in graphene-based system [Eq. (!22|) ]. According to Eq. 
for the maximum of the transition temperature one obtains 

T,"^^" ^ exp{-N)eF = exp(-4)e^ ^ IQ-^f. (45) 

This value is five orders higher than that for the graphene double-layer [Eq. ([2])] and 
by no means does it lead to similar pessimistic conclusions about the feasibility of the 
effect. Of course, in this case the estimate (H5l) may already be quite crude because 
•^A=1) ~ 1/4 is not that small. However, this only indicates and the system could be 
on the verge of applicability of the weak-coupling limit and the actual Tc could be even 
higher than predicted by Eq. (H5l) . 



8. Excitonic condensation in a single- layer graphene 

Besides a double-layer graphene system (Fig. [1]), there were a couple of theoretical 
proposals how excitonic pairing could be realized in a single layer of graphene. We 
briefly discuss them in this section. 

One possibihty [28] would be to apply a strong in-plane magnetic field to a 
single layer of graphene. Such a field acts on the spins of the carriers only, and the 
Zeeman splitting creates electrons with one spin polarization and holes with the opposite 
polarization in an initially neutral sample. As electrons and holes in such a system have 
well-defined Fermi surfaces, the transition temperature can be estimated within the 
large- approximation using Eqs. (fT2|) and (fT3|) of the BCS-like theory. The number of 
species in a single-layer graphene is Ni = NsN^ = 4 and the maximum of the coupling 
constant f[T3|) equals to 

The Zeeman splitting energy ez plays the role of the Fermi energy in such a setup. 
Therefore, for the maximum of the transition temperature one obtains 

T^" ^ exp{-X2To) ez = exp(-8) ~ 3 • IQ-'ez (47) 

The exponential factor exp(— 8) is not as small as exp(— 16) in the double-layer system. 
However, the Zeeman splitting energy cannot be extremely high even for experimentally 
very high magnetic fields B. For S ^ 40 T one can estimate from Eq. (17r|) T^^^ ~ 20mK. 
Alternatively, instead of applying external magnetic field, one could bring the graphene 
sheet in contact with the ferromagnet, in which case the Zeeman splitting would be 
caused by the ferromagnet exchange field [28] . 

In the electrically biased double-layer graphene and a single-layer graphene in 
the parallel magnetic field electrons and holes have a finite density of states at the 
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mutual Fermi level. This way the excitonic pairing is realized though the Cooper 
instability mechanism, and a finite transition temperature formally exists for arbitrary 
small strength of the bare Coulomb interactions (although, as it was demonstrated 
here, even for moderate to strong interactions the transition temperature appears to be 
numerically extremely small). It was conjectured in a number of works [331 EH ES, [36] 
that the requirement of the finite DoS in actually not necessary and that excitonic 
pairing could be possible in an undoped graphene layer, with the chemical potential 
at the Dirac point. The analysis of Refs. [331 EH EHl E6] predicts that opening of the 
excitonic gap at the Dirac point may occur, provided the Coulomb coupling strength 
Ts = e^/{ev) exceeds certain threshold value r*. This value depends on the number of 
fiavors and for = A^^i = 4 equals r* = 2.33, which is quite close to 2.16 in vacuum 



9. Conclusion 

We theoretically studied the possibility of excitonic condensation in a recently proposed 
double-layer graphene system. On the one hand, this system possesses several obvious 
engineering advantages compared to semiconductor double quantum wells, used so far 
experimentally. On the other hand, the properties of the graphene spectrum appear 
to be very unfavorable for realization of excitonic condensation in such a system. 
Additional valley degrees of freedom and suppression of backscattering in graphene 
lead to a much smaller value of the screened Coulomb interaction strength than in a 
system with conventional electron spectrum. It is shown the effect of screening on the 
excitonic pairing is just as strong in the excitonic state as it is in the normal state. For 
this reason, the value of the excitonic gap is found to fully agree with the previously 
obtained [26] estimate for the mean-field transition temperature, the maximum possible 
value ~ 10~''ei? of both being in the ImK range for a perfectly clean system. These 
findings disprove the predictions [22l [23l [27] of the excitonic condensation at room 
temperatures and, unfortunately, render experimental observation of this undoubtedly 
interesting effect in a double-layer graphene system unlikely. 
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